Nonadiabatic corrections to rovibrational levels of H2 
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The leading nonadiabatic corrections to rovibrational levels of a diatomic molecule are expressed in 
terms of three functions of internuclear distance: corrections to the adiabatic potential, the effective 
nuclear mass, and the effective moment of inertia. The resulting radial Schrodinger equation for 
nuclear motion is solved numerically yielding accurate nonadiabatic energies for all rovibrational 
levels of H2 molecule in excellent agreement with previous calculations by Wolniewicz. 
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I. INTRODUCTION 

In the fully nonadiabatic approach the total nonrela- 
tivistic energy of a molecular state is obtained by solv- 
ing the Schrodinger equation with kinetic energy of elec- 
trons and of nuclei on the same footing. This approach 
has been applied to vibrational states of several small di- 
atomic molecules [J 0, H, 0, Q • Much more commonly 
though, the total energy of a molecular state is obtained 
in a two-step procedure based on the Born-Oppenheimer 
(BO) approximation [1,0] in which a separation of elec- 
tronic and nuclear motion is assumed. Namely, in the 
first step, the electronic Schrodinger equation with the 
clamped nuclei Hamiltonian is solved for different nu- 
clear configurations yielding the electronic energy as a 
function of the nuclear coordinates. This function, called 
the potential energy surface (PES), serves as a potential 
for the motion of nuclei in the nuclear Schrodinger equa- 
tion. Solving this equation in the second step yields the 
total rovibronic energy of the system. 

The accuracy of theoretical predictions for molecular 
states, limited due to the BO approximation, can be in- 
creased by including the adiabatic, relativistic, and radia- 
tive corrections without renouncing the notion of PES. 
Further increase in the accuracy requires the nonadia- 
batic effects to be taken into account. A desirable way of 
inclusion of these effects is in terms of a geometry depen- 
dent function, which can be added to PES in the same 
manner as all the other corrections. On one hand, the 
nuclear Schrodinger equation, when solved with such a 
potential, gives the molecular energy levels with spectro- 
scopic precision. On the other, the notion of the PES is 
preserved with all its advantages. 

Several, more or less successful attempts to construct 
such a nonadiabatic correction function for a diatomic 
molecule can be found in literature [1, 0, [H, |Tl|, IT^ . 
Bunker and Moss have derived [sj, in the second or- 
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der of the perturbative expansion, an effective nuclear 
Hamiltonian for the ground electronic state of diatomic 
molecules, in terms of the nonadiabatic potential, vibra- 
tional and rotational masses. Assuming constant vibra- 
tional and rotational masses and neglecting the nona- 
diabatic potential, they obtained values of these effec- 
tive masses for H2 and D2 by fitting to the experimental 
data. Later Schwenke in [l3| used the Bunker and Moss 
[§, effective Hamiltonian to perform ab initio calcu- 
lations of nonadiabatic corrections for H2 and H2O. His 
results for purely vibrational spectrum of H2 differ from 
that of Wolniewicz by about 20% due to the inaccu- 
rate numerical representation of the wave function and, 
what we demonstrate in this work, due to the approxi- 
mate second order nonadiabatic potential of Bunker and 
Moss 0. The accuracy of Schwenke calculations for 
H2O is probably not higher, but clearly demonstrates 
wide applicability of the perturbative approach. Very 
recently Kutzelnigg et al. (ill . [T3 | performed simplified 
calculations of the nonadiabatic correction and both ef- 
fective masses as functions of the internuclear distance in 
HJ and H2. In our recent work [T6l| we have introduced 
nonadiabatic perturbative theory and derived formulae 
for the leading nonadiabatic corrections to energies and 
wave functions. The formula for the nonadiabatic energy 
from that work, although apparently different, is in fact 
equivalent to that of Bunker and Moss 1 9] . Our results for 
rotationless vibrational states have been obtained as the 
expectation value of nonadiabatic corrections with the 
adiabatic wave function. Although numerically accurate, 
due to the neglected third order nonadiabatic corrections 
[see Eq. (1T7|]. our results differed by about 2% from the 
previous calculations by Wolniewicz fT5|| and by Stanke 
et al. [17]. 

In this paper, the nonadiabatic perturbation theory 
has been extended in two directions. Firstly, we general- 
ize the previous derivation to rotational states. Secondly, 
we include the previously missing third order correction, 
which has proved significant. Moreover, we present a 
rigorous formulation of the nonadiabatic perturbative 
theory and include the numerical example of the H2 
molecule. This can be extended to any diatomic molecule 
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and potentially to an arbitrarily large molecule. We de- 
rive formulae valid to all orders, present the leading cor- 
rections of order 0{iJ,~'^), and express them in terms of 
the nonadiabatic correction to the potential and the effec- 
tive i?-dependent nuclear mass and the moment of iner- 
tia. These three functions enter the nuclear Schrodinger 
equation, which can be solved numerically for an arbi- 
trary energy level. As a test of the presented perturba- 
tive theory, we perform calculations of all 301 rovibra- 
tional levels of H2 molecule. We find an excellent 0.1% 
agreement with the accurate nonadiabatic corrections for 
states with the angular momentum J < 10, which were 
obtained by Wolniewicz in [l^, and present for the first 
time results for states with J > 10. 



II. THE ADIABATIC APPROXIMATION 

The total wave function cj) is the solution of the sta- 
tionary Schrodinger equation 



with the Hamiltonian 

H = Hci + Hn , 



(1) 



(2) 



split into the electronic and nuclear parts. In the elec- 
tronic Hamiltonian H^i 



(3) 



In the adiabatic approximation the total wave function 
of an arbitrary molecule 



(7) 



is represented as a product of the electronic wave func- 
tion (pel and the nuclear wave function We note, that 
0oi depends implicitly on the nuclear coordinates R. The 
electronic wave function obeys the clamped nuclei elec- 
tronic Schrodinger equation 



[H,i - £,i{R)] \cj),i) 0, 



(8) 



while the nuclear wave function is a solution to the 
Schrodinger equation in the effective potential generated 
by electrons 



[Hn + £^{R) + £ci{R) - E,] \x) = 0, 



where 



^a(^) = (0cl|ff„|0cl), 



(9) 



(10) 



For the diatomic molecule the nuclear radial equation 
reads 

1 d R^ d 



J(J+ 1) 



+5a(i?)+£cl(i?)-^a 



XjiR) 



(11) 



where J is the rotational quantum number. 



nuclear masses are, by definition, set to infinity, and the 
potential V includes all the Coulomb interactions with 
fixed positions Ra of nuclei. The nuclear Hamiltonian 
involves kinetic energies of all nuclei 



2 Ma 



(4) 



The separation of center of mass motion and the choice of 
the reference frame depends on the particular molecule. 
For example, for a diatomic molecule in the space fixed 
reference frame attached to the geometrical center of two 
nuclei takes the form 



where 



V2 




( ' 


1 




2^n 


\Mb 


Ma 



Vfl • Vol , 



(5) 



(6) 



R = Rab = Ra- Rb, and 1/^,, = 1/Ma + I/Mb is the 
nuclear reduced mass. The last term in Eq. ([5|) vanishes 
for homonuclear diatomic molecules. 



III. PERTURBATIVE FORMALISM 



The total wave function 



(12) 



is the sum of the adiabatic solution and a nonadiabatic 
correction. The nonadiabatic correction S(j)m^ is decom- 
posed into two parts 

(50na = 0cl <5X -I- (5'0na , (13) 

which obey the following orthogonality conditions 



(<^'0na|0ol)cl = 0, 

{6x\x) = 0. 



(14) 

(15) 



The last equation means that normalization of (f) is of the 
form 



1X1'/') = 1. 



The total energy 



(16) 



(17) 
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is the sum of the adiabatic energy i?a and the nonadi- 
abatic correction (5i?na- Using above definitions we pro- 
ceed with the derivation of perturbative formulae. 

The starting point is the Schrodinger equation ([T]) with 
the Hamihonian H, the wave function 0, and the energy 
E decomposed into adiabatic and nonadiabatic parts 

[(i/cl-£ol) + (fcl + i/„-Sa-(5£;„a)]|0cl (x+<5x) +'5>na) = 0. 

(18) 

One rewrites this equation to the form 



{£c\ + Hn — — 5Ena) 

l0ci(x + <5x) + <5V„a) (19) 



and, since ^'(/)na is orthogonal to cjjci, Eq. (|14p . the formal 
solution 



I^Vna) 



1 



+ (£el + ffn - - 6En^.)\6' cf)^^)] , (20) 



is obtained, where the prime in the denominator denotes 
subtraction of the reference state from the Hamiltonian 
inversion. When dx and ^'(/)na on the right hand side are 
neglected, Eq. becomes the leading nonadiabatic 

correction to the wave function. In the next step one 
takes Eq. p8)) and multiplies it from the left by {(f>ci\ 



{(b,i\£ei + Hn-E,-SE^,\(b,i (x + ^x)+'5Vna)ci = 0. (21) 
Since x satisfies Eq. © the above can be simplified to 



i£el + £^. + H,-,-E,)\6x) 

5Ena\X + Sx) - (0cl|i?n|^'0na)Gl 



(22) 



and due to Eq. (fTS]) the solution is 

1 



In the last step, one takes Eq. ([22|) . multiplies it from the 
left by and obtains 



(iJa — £ci — £& — Hn)' 

((0el|i?„|<5'0„a)cl - ^i^nalX + ^X))- (23) 



SEae, = {<Pc\X\Hn\S'(f)nii)- 



(24) 



The set of recursive equations (l20l) . ((23)) . and ((24l) forms 
the perturbative expansion of the wave functions S'cj)^^, 
Sx and energy SEaa- For example, starting from ()24|) one 
gets 



SEn, 



{^clX\Hn 
+ (fcl 



1 



i^nl^cl (X + '^X)) 



i£el - H.iY 

H,,-E,~SE^,)\S'4>n.) 



(25) 



which is the sum of the leading, Eq. ([28)) . and the higher 
order nonadiabatic correction, Eq. (j44p . This perturba- 
tive expansion in general assumes that £gi + — Eg_ 
is small with respect to the electronic excitation energy. 



It is not always true, especially for rovibrational levels 
close to the dissociation threshold. In spite of this fact, 
we claim that each power of £ci+^^n— ^^a in these particu - 
lar matrix elements is at least of the order 0{yj md Hn), 
which we demonstrate in next sections for the leading 
terms (J^^^i^na and 6'^'^^ E-^g. of the nonadiabatic perturba- 
tive expansion. 



A. Second-order nonadiabatic corrections 

In the leading order of perturbative treatment of nona- 
diabatic effects one has 



1 



{£c\ — Hc\) 



(26) 



1 



{Ea — £cl — £a — Hn 
1 



0olX 



7(0ol|-ffn|^'(?!>na)gl , 

(27) 

H„ 0cixV (28) 



The general formula (|28p . following [16|, can be read- 
ily rearranged to a more practical form. From now on we 
consider the homonuclear diatomic two-electron molecule 
and separate out electronic matrix elements from the nu- 
clear ones 



X*x{Hn 0ol 



1 



X* VAX 
V|jX*X 



Hn (t>cl 



(fel - -ffel)' 
1 



Vi 



VAx*vk 



R9el 



i£ci - H.iY 
1 



i£ci - H.iY 

1 



VA'/'ci) (29) 

el 



(fei - H,^Y 



cl 



X*XU{R)-V},[x'x] V\R) 



^kx'KxW'^iR) 



(30) 



where the last equation is the definition of potentials 
U, V" and W*-' . For the S electronic state (pd 

V = n'V, (31) 
W'^ = n' W|| + {5'^ - n' n^) W± , (32) 

where n = R/R, hence 

S^^^Ena = J d^R(^X*XS£na{R) 

+n'n^Wkx*Kxy^\\{R) 

+ iS^^ - n= ) X* K X (R) ]■ (33) 
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The function 



B. Third-order nonadiabatic correction 



SSr^^R) = U{R) + ( I + A ) v(i?) (34) 



is the nonadiabatic correction to the adiabatic energy 
curve £c\{R) + £a{R) and pseudopotcntials U, V, W are: 



U{R) 
V{R) 



1 



1 



1 



1 



(35) 



71 • Vfi0ci ) , (36) 



(fel - i?cl)' 



(37) 



>V^(i?) = — 



1 ((5'J -n^Ti^ ) 



x( V 



R(Pcl 



(fel - i?el)' 



(38) 



In order to simplify the nonadiabatic correction of Eq. 
P3p . one notes that the nuclear wave function x has a 
definite angular momentum 



(39) 



where Yjm are spherical harmonics, thus the nonadia- 
batic correction can be transformed to the form 



(40) 



+X'j 



<5£:na(i?) + A^^^Wj.(i?) 



Let us note, that Eq. can also be expressed in 

terms of an expectation value of an effective nonadia- 
batic Hamiltonian (5iJna 



(41) 



where 



-<5£„a(i?). 



(42) 



Now, the nonadiabatic correction 5x of Eq. ((27)) can be 
conveniently rewritten in terms of i5iJna 



1 



-5H^a\x.j): (43) 



The third order nonadiabatic correction of Eq. ((25|) is 

1 



(fel-i?cl)' " 



0clX 



0olX 



1 



(fol — He 



cbeiSx), (44) 



where Sx is given in Eq. (j27|) . Let us split this sum into 
two parts accordingly 



(45) 



While the first part S'-'^^E'^^^ involves many terms which 
are negligible, since they include the third power of /i„ 
in the denominator, the dominating 0(/i~^) term is 



(fel-i?el)' '^R^^^'^Rx)+0{f^n') (46) 

d3i?n. Vfl(x*x)'5V(i?) + 0(M-'), (47) 



where 



SV{R) = 



1 ^gcl 

2/^2 dR 



(48) 



and this correction is included into S£na of Eq. (|34p . 
which becomes now 

<5£„a(i?) - KiR) + (I + ^) [^(^) + ^^(^)] ■ (49) 

The second term S^'^^'E'^^ can be obtained from the non- 
perturbative solution of the nuclear equation with the 
nonadiabatic Hamiltonian SHna from Eq. (f42|) . Namely, 
for the states with the rotational quantum number J, 
5'^^^E'' takes the form 



S^-^'>Ka - {Xj\SH^a\Sxj) 



(50) 



XJ 



1 



-6Hn 



XJ 



which is interpreted as a second order correction due to 

SHna- 



C. Nuclear equation with effective potentials 



where it is understood, that the derivatives with respect 
to electronic variables of the function x do vanish. 



Instead of calculating the second order correction 
Eq. (|1T|) and the third order correction of Eq. ([50]) . it 
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is more appropriate to solve nonperturbatively the cor- 
responding equation 

[i7„ + £oi{R) + £4R) + SHn^] XJ^E xj , (51) 

where the total energy E is 

E^E, + S''^^En, + S'-''^E^, + ... (52) 

and this is done in this work. Moreover, from comparison 
of subsequent terms of Eq. (|42|) with those of Eq. (fTTjl . 
one concludes that the first term can be interpreted as 
an i?-dependent correction to the nuclear reduced mass 

/in, 



1 



2Ai||(i?) 



^ + W||(i?), 



(53) 



whereas the second term — as an i?-dependent correction 
to the inverse of the moment of inertia (R) 



1 



2fi±iR) 



1 

2 /in 



W±(i?). 



(54) 



With these newly defined functions one can write the 
radial equation as 



1 d R^ 



d 



J(J + 1) 



2i?2 dR fi\\{R) dR 2^^(i?)i?2 
^Ex.j{R). 
where 

y{R) - fel(i?) + fa(i?) + S£UR) 

forms the effective nonadiabatic potential. 



y{R) 



Xj{R) 

(55) 
(56) 



D. 



Asymptotics of the effective masses 



The adiabatic correction £a{R), Eq. (fTU|) . and the 
nonadiabatic correction (5f„a(i?), Eq. do not vanish 
at large internuclear distances. For example, for the large 
atomic separation in the hydrogen molecule, £a.{R) and 
5£na.{R) are equal to me/m-p and —(me/m-pY, respec- 
tively, which corresponds to the first terms in the expan- 
sion of the atomic reduced mass fi — (1/mp -I- l/mo)~^ 
in the electron-nucleus mass ratio 



1 



me /nip 
1 + nie/mp 



rrir.. 



(57) 

Large R asymptotics of the pseudopotentials W|| {R) and 
W±{R) are equal to —mc/rrLp, which is related to the 
change in Eqs. ([55]) and ([5^ of the reduced nuclear mass 
to the reduced mass /xa = (wp + me)/2 of two hydro- 
gen atoms : 



1 



1 



2/i||(oo) 



2fi_i_{oo) 2fiA 
1 1 
m 



m 



me 



= — 1- — + ••• -(58) 



Toe 



E. Evaluation of the wave function derivatives 

The electronic matrix elements in Eq. (j29p involve mul- 
tiple differentiation of the electronic wave function with 
respect to the internuclear distance R, which is difficult 
to calculate directly. Therefore, following Ref. [1^1, we 
rewrite these terms to a more convenient form, where 
differentiation is taken of the Coulomb potential, namely 



1 



{£cl — He 



(59) 



•^^1 = 7? TT^ V^(^)'^oi + 2Vfl(y-£:ei)(60) 



\f-c\ — tle\ ) J 

The derivatives of the potential V 

^ _ I 1 1 1 1 1 
riA riB r2A 7'2S fi2 Rab 

are the following 

1 / riA ^ riB r2A , ?^2B \ R 



(61) 



lA 'IB ' 2A '25 



i?3 



(62) 

V|(y) = 7T [S{nA) + Sine) + S{r2A) + S{f2B)] 

-47r(5(i?), (63) 

and the matrix elements with these operators are readily 
evaluated. The presence of the Dirac delta operators in 
Eq. ([60]) may potentially decrease the accuracy of the 
evaluation of those quantities which contain Vj.0ci- If we 
note that 

Vl{V) = t: [SifiA) + 6{fiB) + Sif2A) + 6if2B)] , (64) 

we can get rid of Dirac deltas by a simple rearrangement 
of the nuclear part of the Hamiltonian to the form 



1 



Hn = -— (Vi^-V,^)-— V^i 



1 



(65) 



The difference in parenthesis collects the terms of 
Eq. ([63|) and (1641) which cancel out, up to the negligi- 
ble —4:Tt5{R) term, so that we can write it down as 



(V| - Vl) = 



(66) 



{£cl — Hcl) 



£cl — Hcl 



with some constant A. In this way the Dirac delta terms 
are eliminated at the expense of employing additional 
basis sets for evaluation of the last resolvent in Eq. (|66p . 
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IV. NUMERICAL CALCULATIONS 

In order to form the radial nonadiabatic equation (|55p . 
apart from the clamped nuclei energy £ei{R) and the adi- 
abatic correction £^{R), it is necessary to know the pseu- 
dopotentials U[R),V{R) + 5V{R),W\\{R), and >V±(i?). 
The evaluation of these functions is the main numeri- 
cal task described here. Calculations were performed at 
80 points including long and very short internuclear dis- 
tances R. At each point several electronic wave func- 
tions of different symmetry were generated. All these 
functions were represented as linear expansions of prop- 
erly symmetrized two-electron basis functions. The basis 
functions were taken in the form of exponentially corre- 
lated Gaussians (EGG) 



Mrur2) = (l + Pi2)(l±i)Sfc 



(67) 



X exp 



where the matrices Afc and vectors Sk contain nonlin- 
ear parameters, 5 per basis function, to be variationally 
optimized. The antisymmetry projector (1 -I- Pi 2) en- 
sures singlet symmetry, the spatial projector (l±i) — the 
gerade or ungerade (— ) symmetry, and the pref- 
actor enforces E states when equal to 1, or 11 states when 
equal to yt — the perpendicular Cartesian component of 
the electron coordinate. 

Six different basis sets have been prepared to repre- 
sent the variety of the electronic wave functions. To en- 
sure high accuracy of the potentials, the basis sets have 
been variationally optimized with respect to pertinent 
goal functions according to the specification in Table HI 



TABLE I: Goal functions used in optimization of the basis 
sets. 



Label Size Symmetry Goal function 



A 


600 




B 


600^" 




C 


1200 




D 


600t 




E 


600 




F 


600 





vj,Xv)-^-—\7Mv) 



' Optimized along with the fixed basis A. 

Particular goal functions have been chosen to reflect 
the contents of the expression the basis set is to be used 
for. The first basis (labeled A), composed of 600 EGG 
functions (|57|) . was employed to expand the X^E+ elec- 
tronic ground state wave function (j^^i- Their nonlinear 



parameters were optimized variationally with respect to 
the clamped nuclei energy with the target accuracy of the 
order of a fraction of microhartree. The bases B, G, and 
D were intended for evaluation of the resolvents present 
in Eqs. ^ and dM])- The two S+ bases (B and D) were 
optimized in the presence of the basis A: the first 600 
terms were taken from wave function and their non- 
linear parameters were kept fixed during the optimiza- 
tion, only the remaining 600 terms were actually opti- 
mized. This ensures that the internal wave function </)ci 
is well represented at every step of optimization. Then, 
the subtraction of the reference state, denoted by the ' 
symbol within the resolvent, was achieved by orthogonal- 
ization of Vu{V) {(pel) to the internal |0ci). In the final 
calculations the three bases A, B, and D were assembled 
together to form a 1800-term basis applied not only 
to evaluate the pertinent resolvents but also to expand 
the external ground state function 0oi. The two ungerade 
bases (E and F) were employed to evaluate the resolvent 
and to form the components of the scalar product in the 
second term of the right hand side of Eq. ([TO]) . 

The adiabatic potential of the nuclear Schrodinger 
equation (|12p was composed of the clamped nuclei energy, 
fci(P), and the adiabatic correction £a{R)- For SeiiR) 
we used the analytic potential constructed by Lach [l^l 
on the basis of the energy points cornputed by Gencek 
from 1200-term EGG wave functions [19| and Sims and 
Hagstrom from Hylleraas wave functions [13]. Their en- 
ergy points where converged up to 13 significant digits. 
The adiabatic correction £a{R) was evaluated as an ex- 
pectation value of the Hamiltonian i/n, Eq. (|10p . 



^a(i?) = 



1 

2^n 



Vcl|(/'cl)cl . 



(68) 



To avoid the cumbersome differentiation of the electronic 
wave function with respect to the internuclear distance 
we replaced the expectation value in the first term on the 
right hand side of Eq. (|68p by an equivalent expression 

(0el|V|| 



which, with the help of Eq 
formed to 



i)ci = -(Vfl^(?!)ci|Vii0oi)ci , (69) 
can be further trans- 



|V^|0. 



l/el 



Imv) 



1 



..of""""),,. 



The formula (|70p . when evaluated with the optimized 
bases A, B, and C, yields the adiabatic correction with 
an accuracy of at least 1 ppm. The adiabatic potential 
curve was then obtained by means of 10-point piecewise 
polynomial interpolation. 

The electronic matrix elements U, V + 5V^ Wn , >V_l en- 
tering Eq. ([55]) were evaluated with the EGG basis sets 
described above, yielding smooth functions of i?. Because 
for the highest vibrational levels the nuclear wave func- 
tions are spread out and the contributions from larger 
internuclear distances are non-negligible, the functions 



7 



IA{R), V(i?), and W(i?) were represented by their asymp- 
totic forms: 



The resulting radial equation 



uo + Uf,/R^ + ^i8/-R^ 

W\\0 + W||i2/i?^^ + W||i4/i?^'', 
W±a + W±12/R^^ + W_L14/i?^**, 



U{R) 
V{R) 
SV{R) 

Wi(i?) « Wio + ^i'±l2/i?'' + ^«±l4/i^'^ (71) 

subject to uq ~ Wj|o = w_lo — — (we/rnp)^ restriction 
(in atomic units). The remaining, free parameters Ui, Vi, 
and Wi were determined by fitting the above functions 
to the calculated points in the range of (6.0, 10.0) bohrs. 
Because at distances i? > 6, the numerical precision of 
the potentials U and V was not high enough, we used 
lower i?-values for the extrapolation. At the origin R = Q 
all the potentials are finite with V ^ R, W|| ^ R^, and 
W_L ~ R^- Numerical results for (5i?na, VV||, and W± are 
shown graphically in Fig. [TJ 




FIG. 1: (Color online) The nonadiabatic potentials (in a.u.): 
(5fna(-R) (bolded, black), W||(7i) (lower dotted, blue), and 
W±{R) (upper dotted, red). The potentials asymptote goes 
at -1/ml = -2.966077 x 10"'^. 



The nonadiabatic correction to energy levels can be 
computed perturbatively from Eq. (j40p as has been done 
in our previous paper 16], or, in a more accurate way, 
the radial equation ([?T|) can be solved for the total nona- 
diabatic energy. This second method is described shortly 
below. The first term including the second order deriva- 
tive is rewritten to the more convenient form 



1 d 



i?2 dR 

_1 A 
~R dR 



152 «15^MT- + W||(i?) — 



2 fin 



2 fin 



OR 



W|(i?) 
R 



(72) 



d 1 d 
dRjAR^dR^^-^^' 



= -2 



E - y{R) - 



W|^(i?) J(J + 1) 



R 



VjiR), (73) 



is solved numerically for the function rjj{R) = Rx.j{R)- 
We used the code developed by W. Johnson, described 
recently in his book [2l|], and modified it to account 
for the dependence of the mass on the internuclear dis- 
tance. In the calculations we used the following con- 
stants '2^]: the proton mass nip ~ 1836.15267247 TOq 
and the energy units conversion factor 1 hartree = 
219474.6313705 cm"!. 



V. RESULTS AND DISCUSSION 

In a molecule, the moving nuclei are 'coated' with elec- 
trons and the amount of additional mass carried by the 
nuclei changes with R. For a homonuclear molecule in 
a purely vibrational state, the effective mass of the nu- 
cleus at given i?, m||(i?), is just twice the reduced mass 
/i||(i?) defined in Eq. l|53)). Analogously, for a rigid ro- 
tating molecule, the effective nuclear mass m±{R) is re- 
lated to fJ.±{R) of Eq. ([54| . Thus, their i?-dependence 
can be determined explicitly from the potentials Wy (R) 
and W±{R), respectively. Fig. [2] illustrates the changes 
in the two effective nuclear masses with the internuclear 
distance in H2. The functions mj|(i?) and m±{R) join 
smoothly the proton mass at the united atom limit with 
the hydrogen atomic mass (mp -I- m^) at the separated 
atoms limit. Interestingly, for R > 2.41 a.u., the effec- 
tive mass m||(i?) is greater than the sum of proton and 
electron masses, reaching rup + 1.6 mo at the maximum 
located near R = 3.8 a.u. 

The radial equation (|73p has been solved for all bound 
states with three versions of the potential y{R): 



yiR) 
yiR) 



Sc\{R), 

£el{R) 

£oi{R) 



£.{R), 

£e.{R)^ 



S£n.{R), 



BO 

adiabatic 
nonadiabatic 



yielding three sets of dissociation energies. The corre- 
sponding dissociation thresholds were —1 hartree in the 
BO approximation, —l + nic/mp hartree in adiabatic ap- 
proximation, and — 1 + rric/mp — (mc/mp)'^ hartree in 
the nonadiabatic level of theory. The results are listed 
in Table |TT1 where for each pair of quantum numbers v 
and J three entries are given (in cm^^): the BO dissoci- 
ation energy, the adiabatic correction, and the nonadia- 
batic correction to the dissociation energy. Thus, the to- 
tal nonrelativistic dissociation energy can be obtained by 
summing up all three entries. The only exception is the 
state with v = 14 and J = 4, for which a nonadiabatic 
level lying just beneath the dissociation threshold has 
been predicted although neither BO nor adiabatic bound 
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states exist. The entry given for this state is the energy 
separation from the nonadiabatic dissociation threshold. 



1838.0 


^.^||(R) 


- 


1837.5 




- 

mp+rrig 


1837.0 






1836.5 


/ / m^(R) 




1836.0 







2 4 6 8 10 

R/a.U. 



FIG. 2: (Color online) The ii-dependence of the effective nu- 
clear masses Tn||(iZ) and m±{R) (in a.u.). 

Our nonadiabatic corrections agree very well with 
those computed by Wolniewicz for rovibrational 
states of J < 10. For all rotational states of the lowest 
vibrational level the difference is merely 0.0002 cm~^ or 
less. In all the cases the difference is smaller than 0.1%, 
which corresponds to the magnitude of the neglected 
higher order terms of relative order C'(l//Xn). This is 
the first confirmation of the Wolniewicz's results for 
J > rovibrational states and, simultaneously, a numer- 
ical validation of the nonadiabatic perturbation theory 
presented in this work. Let us note, that the nonadia- 
batic corrections to the dissociation energy from Table HIl 
differ in sign from the corrections to the adiabatic energy 
of [l^ and their absolute values differ by the constant 
l/Wp = 0.06509787 cm~^ corresponding to the asymp- 
totic value of the nonadiabatic potential SSnniR), dis- 
cussed in Subsec. IIIIDI 

VI. SUMMARY 

We have presented the nonadiabatic perturbation the- 
ory applicable to any molecule in an arbitrary rovibra- 
tional state. The leading nonadiabatic corrections for 
the diatomic molecule are expressed in terms of three 
i?-dependent functions depicted in Fig. [TJ the nuclear 
reduced mass in Eq. (j53p . the moment of inertia in 
Eq. ([5^ . and the correction Eq. to the adiabatic 
potential, which enter the radial Schrodinger equation 
for nuclei Eq. ((5T|) . This equation can be solved for an 
arbitrary molecular states by standard numerical meth- 



ods [2l|. Although representation of the nonadiabatic 
correction by pseudopotentials has been studied previ- 
ously (see and references therein), we have rig- 
orously derived new expressions for W|| , W±_ , and SSna 
functions of the internuclear distance, which give nona- 
diabatic corrections with 0{mc/ fj.n) accuracy. Moreover, 
we have performed explicit numerical calculations for the 
simplest example of II2 molecule in order to verify the ob- 
tained perturbative formulae. All the electronic matrix 
elements with differentiation of the wave function over 
the internuclear distance were rewritten in a convenient 
form involving differentiation of the Coulomb potential. 
Such an approach enables achieving much higher numer- 
ical precision even for the well known adiabatic correc- 
tion. The final accuracy of all three nonadiabatic func- 
tions is limited only by the neglected higher order terms, 
namely 0(TOc//Xn) relative to the leading order, and re- 
sults for rovibrational levels agree within this uncertainty 
with former results of Wolniewicz [l5] . Having accurate 
nonrclativistic energies one can include relativistic and 
QED corrections, which become significant for the states 
close to the dissociation threshold. Relativistic correc- 
tions for the electronic ground state of H2 are known 
accurately for a wide range of internuclear distances [1^ . 
Their large R asymptotics, including QED corrections, is 
presently investigated by the Jeziorski group [23| . Com- 
bining all the knowledge together would enable achieving 
at least an order of magnitude increase in the precision 
of theoretical prediction for all molecular states of II2 . 

Analogous calculations can be performed for D2 and 
HD molecules. It is worth noting, that in the former case 
there is no need to recompute the pseudopotentials and 
only a proper rescaling due to different reduced masses 
is required. The calculations on the latter system, due 
to difference in mass of the proton and deuteron, would 
involve additional correction to SS^a. coming from the 
last term of Eq. ([5]). While this perturbative approach 
can be further extended to larger diatomic molecules, it 
would be more challenging to investigate three-atomic 
molecules such as Hj^ — a system of great astrophysical 
interest. We expect no principal difficulties in such cal- 
culations except for much more increased demands in 
computer resources needed to perform optimization of 
pertinent wave functions. 
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TABLE II: The dissociation energy of the rovibrational states of H2 (in cm~^). For each pair of vibrational (v) and rotational 

(J) quantum numbers, three entries are given: the BO dissociation energy, the adiabatic correction, and the nonadiabatic 
correction. The sum of the three numbers gives the total nonrelativistic dissociation energy of the {v, J) state. 



v\J 





1 


2 


3 


4 


5 


6 


7 





36112.5927 


35994.0372 


35758.0143 


35406.6660 


34943.1216 


34371.4040 


33696.3152 


32923.3094 





5.7711 


5.8348 


5.9611 


6.1481 


6.3931 


6.6924 


7.0418 


7.4367 





0.4339 


0.4406 


0.4539 


0.4739 


0.5005 


0.5338 


0.5736 


0.6200 


1 


31949.1892 


31836.5534 


31612.3337 


31278.6001 


30838.3755 


30295.5437 


29654.7385 


28921.2187 


1 


7.1740 


7.2311 


7.3443 


7.5119 


7.7310 


7.9983 


8.3096 


8.6605 


1 


1.2704 


1.2761 


1.2874 


1.3044 


1.3271 


1.3555 


1.3896 


1.4294 


2 


28021.4345 


27914.5889 


27701.9154 


27385.4159 


26968.0130 


26453.4605 


25846.2343 


25151.4125 


2 


8.3336 


8.3840 


8.4839 


8.6314 


8.8240 


9.0583 


9.3304 


9.6358 


2 


2.0271 


2.0318 


2.0414 


2.0557 


2.0749 


2.0989 


2.1279 


2.1618 


3 


24324.4498 


24223.3066 


24022.0055 


23722.4851 


23327.5744 


22840.9061 


22266.8106 


21610.1998 


3 


9.2420 


9.2855 


9.3714 


9.4982 


9.6632 


9.8633 


10.0945 


10.3525 


3 


2.7088 


2.7127 


2.7206 


2.7326 


2.7485 


2.7686 


2.7929 


2.8214 


4 


20855.2072 


20759.7245 


20569.7150 


20287.0583 


19914.4976 


19455.5548 


18914.4286 


18295.8815 


4 


9.8890 


9.9251 


9.9966 


10.1017 


10.2378 


10.4019 


10.5902 


10.7985 


4 


3.3187 


3.3219 


3.3283 


3.3379 


3.3508 


3.3670 


3.3868 


3.4100 


5 


17612.7145 


17522.9036 


17344.2107 


17078.4626 


16728.3253 


16297.2225 


15789.2368 


15209.0018 


5 


10.2629 


10.2915 


10.3476 


10.4298 


10.5356 


10.6618 


10.8049 


10.9606 


5 


3.8576 


3.8599 


3.8647 


3.8719 


3.8816 


3.8938 


3.9087 


3.9262 


6 


14598.2891 


14514.2237 


14346.9982 


14098.3935 


13771.0098 


13368.1874 


12893.9129 


12352.7160 


6 


10.3525 


10.3729 


10.4130 


10.4712 


10.5449 


10.6313 


10.7268 


10.8271 


6 


4.3218 


4.3232 


4.3262 


4.3306 


4.3365 


4.3440 


4.3531 


4.3639 


7 


11815.9479 


11737.7779 


11582.3242 


11351.3310 


11047.3466 


10673.6480 


10234.1517 


9733.3170 


7 


10.1486 


10.1606 


10.1839 


10.2169 


10.2573 


10.3020 


10.3476 


10.3901 


7 


4.7029 


4.7031 


4.7037 


4.7044 


4.7055 


4.7067 


4.7081 


4.7095 


8 


9272.9561 


9200.9278 


9057.7445 


8845.1263 


8565.5890 


8222.3732 


7819.3619 


7360.9906 


8 


9.6483 


9.6517 


9.6577 


9.6648 


9.6707 


9.6726 


9.6670 


9.6500 


8 


4.9838 


4.9824 


4.9796 


4.9752 


4.9692 


4.9614 


4.9516 


4.9393 


9 


6980.5984 


6915.0837 


6784.9238 


6591.8323 


6338.3203 


6027.6323 


5663.6708 


5250.9187 


9 


8.8590 


8.8537 


8.8422 


8.8232 


8.7945 


8.7533 


8.6963 


8.6196 


9 


5.1364 


5.1325 


5.1246 


5.1124 


5.0958 


5.0743 


5.0472 


5.0138 


10 


4955.2699 


4896.8118 


4780.7739 


4608.8912 


4383.7137 


4108.5514 


3787.4126 


3424.9438 


10 


7.8021 


7.7882 


7.7597 


7.7151 


7.6525 


7.5691 


7.4616 


7.3260 


10 


5.1160 


5.1082 


5.0925 


5.0684 


5.0354 


4.9926 


4.9386 


4.8720 


11 


3220.0418 


3169.4253 


3069.1017 


2920.8723 


2727.3991 


2492.1652 


2219.4373 


1914.2414 


11 


6.5140 


6.4918 


6.4469 


6.3777 


6.2821 


6.1573 


5.9996 


5.8045 


11 


4.8566 


4.8429 


4.8151 


4.7725 


4.7140 


4.6379 


4.5418 


4.4224 


12 


1806.9489 


1765.3260 


1683.0605 


1562.1049 


1405.3733 


1216.7368 


1001.0464 


764.2101 


12 


5.0372 


5.0068 


4.9451 


4.8505 


4.7203 


4.5506 


4.3360 


4.0685 


12 


4.2657 


4.2429 


4.1965 


4.1254 


4.0271 


3.8984 


3.7345 


3.5281 


13 


760.3903 


729.5279 


668.9437 


580.9383 


469.0254 


338.0600 


194.5461 


47.4825 


13 


3.3933 


3.3526 


3.2697 


3.1417 


2.9631 


2.7253 


2.4123 


1.9886 


13 


3.2221 


3.1850 


3.1095 


2.9927 


2.8294 


2.6113 


2.3232 


1.9317 


14 


141.7951 


124.7523 


92.3077 


48.0033 










14 


1.5343 


1.4739 


1.3479 


1.1416 










14 


1.5847 


1.5226 


1.3933 


1.1825 


0.0887+ 









' This state appears as a resonance in BO and adiabatic approximations. The entry is a dissociation energy of this 
nonadiabatic level. 
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TABLE II: continued 



v\J 


8 


9 


10 


11 


12 


13 


14 


15 





32058.3583 


31107.8190 


30078.3065 


28976.5781 


27809.4302 


26583.6105 


25305.7462 


23982.2860 





7.8718 


8.3418 


8.8411 


9.3642 


9.9055 


10.4596 


11.0212 


11.5852 





0.6728 


0.7322 


0.7981 


0.8704 


0.9491 


1.0343 


1.1260 


1.2240 


1 


28100.7394 


27199.4234 


26223.6397 


25179.8929 


24074.7260 


22914.6378 


21706.0160 


20455.0847 


1 


9.0458 


9.4605 


9.8991 


10.3562 


10.8264 


11.3045 


11.7850 


12.2629 


1 


1.4749 


1.5263 


1.5834 


1.6464 


1.7154 


1.7903 


1.8713 


1.9584 


2 


24374.5496 


23521.5539 


22598.5709 


21611.8782 


20567.7938 


19472.5994 


18332.4783 


17153.4682 


2 


9.9697 


10.3271 


10.7026 


11.0910 


11.4869 


11.8851 


12.2802 


12.6671 


2 


2.2008 


2.2449 


2.2942 


2.3488 


2.4089 


2.4744 


2.5456 


2.6225 


3 


20876.4463 


20071.2651 


19200.6029 


18270.5391 


17287.1990 


16256.6838 


15185.0135 


14078.0862 


3 


10.6327 


10.9300 


11.2393 


11.5554 


11.8729 


12.1865 


12.4908 


12.7804 


3 


2.8543 


2.8917 


2.9337 


2.9804 


3.0320 


3.0885 


3.1500 


3.2167 


4 


17605.1247 


16847.7045 


16029.3982 


15156.1198 


14233.8415 


13268.5282 


12266.0890 


11232.3425 


4 


11.0221 


11.2561 


11.4954 


11.7348 


11.9690 


12.1923 


12.3995 


12.5849 


4 


3.4370 


3.4677 


3.5023 


3.5409 


3.5836 


3.6304 


3.6814 


3.7365 


5 


14561.5917 


13852.4141 


13087.1119 


12271.4764 


11411.3755 


10512.6972 


9581.3093 


8623.0355 


5 


11.1244 


11.2915 


11.4566 


11.6146 


11.7601 


11.8873 


11.9906 


12.0642 


5 


3.9466 


3.9699 


3.9962 


4.0254 


4.0576 


4.0926 


4.1304 


4.1704 


6 


11749.5653 


11089.7678 


10378.8781 


9622.6208 


8826.8275 


7997.3920 


7140.2429 


6261.3368 


6 


10.9279 


11.0242 


11.1110 


11.1829 


11.2344 


11.2597 


11.2528 


11.2071 


6 


4.3763 


4.3903 


4.4058 


4.4226 


4.4406 


4.4592 


4.4779 


4.4957 


7 


9176.0486 


8567.6048 


7913.5170 


7219.5226 


6491.5173 


5735.5270 


4957.7045 


4164.3543 


7 


10.4250 


10.4474 


10.4524 


10.4343 


10.3876 


10.3060 


10.1831 


10.0114 


7 


4.7108 


4.7118 


4.7119 


4.7108 


4.7076 


4.7012 


4.6902 


4.6724 


8 


6852.1602 


6298.1569 


5704.5851 


5077.3199 


4422.4828 


3746.4465 


3055.8806 


2357.8507 


8 


9.6172 


9.5638 


9.4847 


9.3741 


9.2260 


9.0337 


8.7892 


8.4833 


8 


4.9241 


4.9053 


4.8820 


4.8528 


4.8160 


4.7692 


4.7092 


4.6313 


9 


4794.3657 


4299.4465 


3771.9991 


3218.2507 


2644.8437 


2058.9184 


1468.2878 


881.7751 


9 


8.5188 


8.3891 


8.2250 


8.0203 


7.7678 


7.4587 


7.0814 


6.6193 


9 


4.9730 


4.9232 


4.8624 


4.7880 


4.6962 


4.5819 


4.4375 


4.2515 


10 


3026.3806 


2597.5217 


2144.7391 


1675.0460 


1196.2612 


717.3517 


249.1650 




10 


7.1577 


6.9514 


6.7003 


6.3963 


6.0280 


5.5786 


5.0193 




10 


4.7905 


4.6909 


4.5694 


4.4200 


4.2343 


3.9995 


3.6930 




11 


1582.3694 


1230.4417 


866.0734 


498.2452 


138.1669 








11 


5.5658 


5.2756 


4.9227 


4.4892 


3.9416 








11 


4.2752 


4.0939 


3.8689 


3.5852 


3.2148 








12 


513.3908 


257.4738 


8.3253 












12 


3.7358 


3.3166 


2.7633 












12 


3.2682 


2.9354 


2.4878 
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TABLE II: continued 



v\J 


16 


17 


18 


19 


20 


21 


22 


23 





22619.4577 


21223.2371 


19799.3296 


18353.1603 


16889.8728 


15414.3355 


13931.1527 


12444.6818 





12.1465 


12.7005 


13.2425 


13.7681 


14.2729 


14.7525 


15.2029 


15.6196 





1.3286 


1.4396 


1.5572 


1.6814 


1.8122 


1.9497 


2.0939 


2.2451 


1 


19167.8656 


17850.1517 


16507.4922 


15145.1868 


13768.2885 


12381.6135 


10989.7585 


9597.1235 


1 


12.7332 


13.1911 


13.6317 


14.0502 


14.4420 


14.8023 


15.1262 


15.4086 


1 


2.0516 


2.1512 


2.2570 


2.3693 


2.4881 


2.6134 


2.7453 


2.8838 


2 


15941.4274 


14702.0131 


13440.6714 


12162.6372 


10872.9422 


9576.4337 


8277.7999 


6981.6067 


2 


13.0406 


13.3957 


13.7272 


14.0300 


14.2988 


14.5281 


14.7122 


14.8446 


2 


2.7052 


2.7938 


2.8884 


2.9890 


3.0955 


3.2080 


3.3263 


3.4502 


3 


12941.6497 


11781.2863 


10602.4104 


9410.2757 


8209.9952 


7006.5729 


5804.9493 


4610.0657 


3 


13.0501 


13.2944 


13.5077 


13.6843 


13.8181 


13.9026 


13.9307 


13.8939 


3 


3.2886 


3.3657 


3.4479 


3.5352 


3.6272 


3.7237 


3.8240 


3.9274 


4 


10172.9970 


9093.6453 


7999.7724 


6896.7794 


5790.0229 


4684.8764 


3586.8202 


2501.5762 


4 


12.7427 


12.8672 


12.9520 


12.9905 


12.9756 


12.8991 


12.7513 


12.5207 


4 


3.7957 


3.8587 


3.9253 


3.9948 


4.0665 


4.1392 


4.2113 


4.2803 


5 


7643.6478 


6648.8775 


5644.4446 


4636.1112 


3629.7676 


2631.5640 


1648.1214 


686.8842 


5 


12.1016 


12.0965 


12.0417 


11.9291 


11.7496 


11.4918 


11.1414 


10.6779 


5 


4.2124 


4.2556 


4.2990 


4.3412 


4.3801 


4.4126 


4.4344 


4.4383 


6 


5366.6735 


4462.3376 


3554.5745 


2649.9163 


1755.3887 


878.8631 


29.7112 




6 


11.1158 


10.9712 


10.7645 


10.4851 


10.1192 


9.6476 


9.0402 




6 


4.5115 


4.5235 


4.5291 


4.5249 


4.5055 


4.4629 


4.3836 




7 


3361.9951 


2557.4749 


1758.1694 


972.3301 


209.7400 








7 


9.7827 


9.4869 


9.1113 


8.6376 


8.0373 








7 


4.6449 


4.6038 


4.5430 


4.4532 


4.3193 








8 


1660.0067 


970.9237 


300.7647 












8 


8.1037 


7.6332 


7.0443 












8 


4.5286 


4.3911 


4.2015 













9 
9 
9 



309.8952 
6.0447 
4.0047 



13 



TABLE II: continued 



v\J 


24 


25 


26 


27 


28 


29 


30 


31 





10959.0542 


9478.2014 


8005.8846 


6545.7303 


5101.2727 


3676.0048 


2273.4444 


897.2218 





15.9981 


16.3340 


16.6221 


16.8572 


17.0329 


17.1423 


17.1769 


17.1261 





2.4031 


2.5682 


2.7404 


2.9199 


3.1066 


3.3006 


3.5019 


3.7105 


1 


8207.9415 


6826.3161 


5456.2676 


4101.7929 


2766.9437 


1455.9363 


173.3131 




1 


15.6440 


15.8264 


15.9492 


16.0045 


15.9831 


15.8734 


15.6600 




1 


3.0290 


3.1807 


3.3389 


3.5032 


3.6734 


3.8488 


4.0283 




2 


5692.3452 


4414.4956 


3152.6140 


1911.4560 


696.1630 








2 


14.9184 


14.9254 


14.8559 


14.6980 


14.4356 








2 


3.5795 


3.7135 


3.8516 


3.9923 


4.1340 








3 


3426.9550 


2260.8735 


1117.5022 












3 


13.7827 


13.5850 


13.2850 












3 


4.0327 


4.1379 


4.2403 












4 


1435.3165 


395.0112 














4 


12.1920 


11.7435 














4 


4.3427 


4.3926 















